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?— ( Abstract. We present a numerical study of the effect of inertia 



< 



on suspended particles moving through a geometric constriction 
to understand the behavior of particles in analogous microfluidic 

■^ settings, such as pinched flow fractionation devices. The parti- 

CN cles are driven by a constant force in a quiescent fluid, and the 

^^ constriction (the pinching gap) corresponds to the space between 

r^ a plane wall and a second, fixed sphere of the same size (the ob- 

J>-j stacle). The results show that, due to the effect of inertia and 

''O geometric constriction (individual or combined) the particles at- 

^ tain smaller surface-to-surface separations from the obstacle. We 

rr^ . then relate the minimum surface-to-surface separation to the ef- 

• feet that the short-range repulsive non-hydrodynamic interactions 

Q (such as solid-solid contact due to surface roughness, electrostatic 

'^ double layer repulsion, etc.) would have on the particle trajec- 

^~, tories. In particular, using a simple hard-core repulsive potential 

i-G model for such interactions, we infer that the particles would ex- 

I '~ < perience larger lateral displacements moving through the pinching 

gap, as inertia increases and/or the aperture created by the con- 

^-H striction decreases. Thus, separation of particles based on the dif- 

^ ferences in density is in principle possible, owing to differences in 

)s^ inertia associated with them. However, we observe in the simul- 

p — taneous presence of pinching wall and significant inertia, that the 

\^ wall reduces inertia effects and hinders separation. 

o 
en 

> 1. Introduction 

/S The deterministic motion of spherical particles through geometric 

d constrictions is at the core of a number of microfluidic unit opera- 

tions and devices. Some prominent examples include separation sys- 
tems, e.g., based on pinched flow fractionation, [1] particle focusing 
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2 A SPHERE THROUGH A CONSTRICTION 

methods, [2] such as geometrical focusing, p] and constricted microchan- 
nels and micromodels of porous media used to study particle deposition 
and clogging. [H [5] In all cases, the interactions of the suspended parti- 
cles with the confining walls of a microfluidic channel play a crucial role 
in determining their motion and fate. The nature of the interactions 
can be hydrodynamic, particularly important when particle size is com- 
parable to the confinement, or non-hydrodynamic, such as solid-solid 
contact due to surface roughness. 

The motion of suspended particles in microfluidic systems is in many 
cases described based on the fluid flow streamlines obtained in the 
absence of particles. Rigorously, such an explanation is valid only 
for point particles, because it neglects the effect of particle size and 
particle-wall hydrodynamic interactions. Moreover, the streamlines are 
calculated in the Stokes regime, neglecting fluid inertia. Particle inertia 
is also ignored and the effect of non-hydrodynamic interactions, such 
as solid-solid contact, is typically incorporated indirectly. For example, 
Jain and Posner model the effect of roughness on particle dispersion 
in order to estimate the relevant statistical properties of the separa- 
tion process, such as its resolution. [6J However, there is ample evidence 
that hydrodynamic interactions, inertia effects, and non-hydrodynamic 
repulsive forces significantly affect the individual particle trajectories 
through narrow constrictions, leading to the separation process. [3 [HI [8] 

Consider the case of pinched flow fractionation (PFF), in which 
the separative lateral displacement experienced by particles of unequal 
size is achieved by passing the suspension through a constriction, fol- 
lowed by a sudden expansion. [U [TUl HU El H^l US] Although the orig- 
inal explanation offered for the separative lateral displacement was 
based on particle motion superimposed on the streamlines in Stokes 
fiow,[T] particle size and hydrodynamic interactions cannot be ignored 
in PFF due to the small particle-wall separations that take place near 
the constriction. [7f 18] It has been also shown recently, that inertia 
effects (manifesting themselves as 'lift' forces) might play a crucial 
role in PFF.p] In some particle focusing methods inertial effects are 
also present, and even exploited as a suspension of particles moves 
through a series of constrictions. [HI |2] In fact, the asymmetric trajec- 
tories observed in focusing devices in which the constriction is fore- 
aft symmetric[3] indicates the presence of irreversible forces, i.e., in- 
ertial forces or non-hydrodynamic interactions between the particles 
and the channel wall. Otherwise, the trajectories should be sym- 
metric and reversible, consistent with Stokes equations. In the case 
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of clogging in microchannels, it has been argued that particle depo- 
sition involves the interplay between hydrodynamic as well as non- 
hydrodynamic interactions. |1] The hydrodynamic interaction of parti- 
cles with constrictions/channel walls leads to particle trajectories that 
deviate from flow streamlines, and if a particle reaches a critical sepa- 
ration from the wall, non-hydrodynamic attractive interactions result 
in particle deposition. [U [3] 

In this work, we investigate the motion of a spherical particle pass- 
ing through a constriction created by a spherical obstacle of the same 
size and a plane wall (see figure Q . We are interested in the minimum 
surface-to-surface separation occuring between the particle and the ob- 
stacle during the motion, since it is the relevant length scale to assess 
the relative importance of short-range non-hydrodynamic interactions 
between them. Specifically, we show that the minimum separation 
decreases as the magnitude of inertia increases and/or the aperture 
created by the constriction decreases. The inferences drawn from these 
results are consistent with the observations reported in previous ex- 
periments, in that, the lateral displacement of a given size of particles 
increases as the aperture of the constriction becomes smaller. [15] In 
addition, the computational nature of the work presented here enables 
us to isolate the effects of fluid and particle inertia by adjusting both 
the particle-to-fiuid density ratio as well as the force acting on the par- 
ticle. In comparison, the ratio of the densities was moderate in the 
experiments reported in Ref. [T3] (0(10°) — O(IO^)), and hence the 
effects of particle and fluid inertia were coupled. 

The article is organized as follows: in the following section, we define 
the geometry of the system under investigation and describe the param- 
eter space relevant to the problem in terms of dimensionless numbers. 
In ^ we introduce a hard-core model for non-hydrodynamic repulsive 
interactions, and discuss its effect on particle trajectories in the ab- 



sence (3.1) and presence (3.2) of inertia. In q4l we present the results 



of our simulations in the presence of a pinching wall only (4.1), in the 



presence of inertia only (4.2), and both in the presence of a pinching 
wall as well as significant inertia (4.2.1). Finally in ^we summarize 



the inferences drawn from this work. 



2. System definition and parameter space 

A schematic view of the system under consideration is shown in 
figure [TJ A spherical particle suspended in a quiescent fluid moves past 
a spherical obstacle of the same radius a (diameter d) under the action 
of a constant force. The x-axis is set in the direction of the constant 
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Figure 1. (a) System, simulation box and boundary 
conditions. The schematic is to scale, i.e., the box is 16 
obstacle/particle radii in length. The top wall creates the 
pinching gap. In the absence of walls, we use the periodic 
boundaries for the box, and the obstacle is at the center 
of the box. (b) The part of the simulation box indicated 
with dashed lines in (a) is enlarged here. The initial 
offset bin, the final offset bout and the aperture of the 
pinching gap D are indicated. The dashed circle depicts 
the position of the particle when the particle-obstacle 
separation is minimum (in the absence of inertia), (c) 
Enlarged view of the particle (dashed) and the obstacle 
(solid) surfaces showing the minimum separation a^min 
between them. 



(C) 



force, and the far upstream (downstream) - i.e., x — )■ — oo {x — )■ +oo) 
- position of the particle along the y-axis is the initial offset bin {final 
offset bout)- A plane wall perpendicular to the positive y-axis creates a 
wall-obstacle pinching gap of minimum aperture D. The motion of a 
particle initially located in the xy-plane is confined to that plane, due to 
the symmetry of the problem. The minimum value of the dimensionless 
surface-to-surface separation between the particle and the obstacle for 
a given trajectory is ^min (the corresponding dimensional separation is 
o-Cmin as shown in figure [I]^c)). 

Numerical simulations are performed using the lattice Boltzmann 
method {susp3d}.^El [TTJ [TH] The boundary conditions imposed in the 
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simulations are shown in figure U[a): no-slip at all fluid-solid inter- 
faces, periodic boundary conditions in x and z directions, as well as in 
^/-direction in the absence of the pinching wall. The size of the simu- 
lation box is 16a in all directions (particle radius a = 5 lattice units). 
We vary the aperture of the pinching gap D by translating the obsta- 
cle center along the y direction. The initial position of the particle, 
upstream to the obstacle (xq ~ —7a), is given by the initial offset bin, 
and the final offset is calculated at the symmetric position with respect 
to the y-axis. Here we note that the code implementing the lattice 
Boltzmann method is designed to construct a channel, such that two 
walls are simultaneously added on the opposite faces of the simulation 
box (for example, see the top and bottom walls in figure [l]( a)). How- 
ever, preliminary simulation tests showed that in a sufficiently large 
box (i.e., 16a in x,y and z directions in this work), the effect of the 
wall opposite to the pinching wall is negligible. 

The parameter space for our investigation is characterized by three 
dimensionless numbers: the Reynolds number Re, quantifying fluid 
inertia, the Stokes number St, quantifying particle inertia, and the 
ratio of the particle diameter to the aperture of the pinching gap d/D, 
quantifying the extent of pinching. We define the Reynolds and Stokes 
numbers based on the bulk (Stokes) settling velocity and the radius of 
the particle, 

(1) 

(2) St 

where Tp = m/67cpa is the inertial relaxation time for an isolated sus- 
pended particle, [19j tj = a /U stokes is the characteristic time of the flow 
around the particle, v is the kinematic viscosity of the fluid {v = 1/6 
in lattice units throughout this worl|j, g is the acceleration due to 
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It has been reported that using kinematic viseosity z/ = 1/6 (in lattice units) has- 
tens the simulations as well as reduces numerical errors [Ladd, A. J. C. and Verberg, 
R. (2001), J Stat Phys, Vol.104, Nos. 5/6, pp. 1191-1251], providing a basis for our 
choice. Also, we note that the accuracy of the lattice-Boltzmann method based on 
the inherent compressibility, is reported to be 0{Ma^) in the literature [e.g., Ladd, 
A. J. C. and Verberg, R. (2001), J Stat Phys, Vol.104, Nos. 5/6, pp. 1191-1251 
and Hazi, G. (2003), Phys Rev E, Vol. 67, article 056705], where Ma is the Mach 
number associated with the system. For our system, the largest instantaneous par- 
ticle velocities observed in the simulations are O(10~^ lattice units). Given that 
the LB-code is D3Q19, velocity of sound is given by c^ = 1/3 in lattice units. Thus, 
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the constant external force acting on the particlen and pp, p/ are the 
densities of the particle and the fluid, respectively. We cover a range 
of Re and St values from negligible (O(10~'^)-O(10~'^)) to appreciable 
(O(O.l)-O(l)) by varying Pp/pf and g, and use four values of the d/D 
ratio: 0.588, 0.5, 0.4 and 0.333. 

In general, the trajectories followed by the particles for different ini- 
tial offsets bin, and thus the minimum dimensionless separations ^mm, 
are functions of the three dimensionless numbers: Re, St, d/D. We 
estimate the effect of short-range non-hydrodynamic interactions on a 
particle trajectory by comparing the range of the interactions with the 
corresponding minimum separation attained by the particle along that 
trajectory. Further, the particle trajectories in the Stokes regime are 
symmetric about the y-axis in the absence of non-hydrodynamic inter- 
actions. Therefore, the symmetry (or equivalently, asymmetry) exhib- 
ited by a particle trajectory also serves a diagnostic purpose: any asym- 
metry present in the trajectories obtained in our simulations is indica- 
tive of deviations from the Stokes regime (i.e. , significant inertia) O [22] • 
The presence of such asymmetry can be determined by comparing the 
initial and final offsets. Thus, in the next sections, we use the minimum 
separation (^mm) as well as the offsets (6j„ and bout) as representative 
variables to investigate the effect of non-hydrodynamic interactions, 
in the presence of inertia and pinching on the particle trajectories. 
Throughout this article, for brevity and clarity, we present results cor- 
responding to bin = 2a. In the supplementary information, we present 
our results for the entire range of 6j„ considered (1.6a-2.2a), and show- 
ing the minimum separation {C,min) as a function of the offsets 6j„ and 
bout for all cases in terms of the aperture of the pinching gap and inertia. 

3. Hard-core model for non-hydrodynamic interactions 

3.1. Stokes regime and the critical offset. We have shown in pre- 
vious work that the minimum separation ^^m decreases dramatically 
with the initial offset bin, especially when bin becomes comparable to 
the obstacle radius (or, the mean of particle and obstacle radii for un- 
equal particle and obstacle sizes). j221 [21] An initial offset of 6j„ ~ a. 



Ma ^ 0(10 "■) and errors related to inherent compressibility of the method are 
0(10-4). 

We note that, in our simulations we specify the constant external acceleration 
acting on the particle, which is multiplied by the particle mass to compute the 
net force acting on it. Therefore, the expression for the Stokes settling velocity in 
equation [T] contains just the particle density Pp, and not the difference between the 
particle and fluid densities, Ap = pp — pf. 



A SPHERE THROUGH A CONSTRICTION 



K < K 


^_^_^ 




^o„.= ^» 


K = K __--/" 


^».= ^. 


/ 


\ \ 
\ 1 


^„„.= ^. 




Obstacle^\ ^^-~ 
\ 
\ 

\ 


1 

/ 1 
/ 1 

/ 

/ 
/ 
y 





Figure 2. Trajectories corresponding to the three pos- 
sible types of particle-obstacle 'collisions', adapted from 
Ref [H]. The particle travels from left (incoming half of 
the trajectory) to right (outgoing half of the trajectory). 
The dashed circle represents the excluded volume inac- 
cessible to the particle centers. 



for example, results in a minimum separation ^min ~ O(10~^). There- 
fore, for a particle with its radius in the millimeter range, the cor- 
responding minimum separation amounts to a few microns, while for 
micrometer-sized particles, the predicted minimum separations are at 
atomic length-scales. In reality, micrometer-sized particles exhibit vari- 
ous repulsive non-hydrodynamic interactions at separations 0(100 nm), 
such as, electrostatic double layer repulsion, [25j solid-solid contact due 
to surface roughness [261 127] and steric repulsion. [28] Following previous 
studies, we use a simple model that describes these interactions with a 
hard-core repulsive potential. [221 EHl EHl EU |32] In this approximation, 
if the dimensionless range of the repulsive interaction is e, any surface- 
to-surface separation ^ < e is not possible, since it is inside the core of 
the interaction. Consequently, if a particle trajectory attains a mini- 
mum separation E,min in the absence of non-hydrodynamic interactions, 
by comparing it with e we investigate the effect of these interactions. 
Moreover, we assume that the hard-core repulsion is the only effect 
of the potential, in that, the hydrodynamic mobility of the suspended 
particle is not affected by the interaction. [221 ESI |30l |3T1 [32] 

In the absence of inertia, there is a one-to-one relationship between 
the minimum separation and the initial offset. [221 [23] We can thus de- 
fine a critical offset he-, as the final offset resulting from a minimum sep- 
aration equal to the range of the repulsive interaction i.e., ^^m = e (see 
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Figure 3. In presence of significant inertia (in this par- 
ticular case, particle inertia, St = 1.11), the trajectories 
are asymmetric even in the absence of non-hydrodynamic 
interactions. Such asymmetry leads to two critical off- 
sets, viz.- bcAn, hrout as shown. 
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figure^. In previous work, we have established that the critical off- 
set he completely characterizes the asymmetry caused by the repulsive 
interaction in the particle trajectory [221 E]: any trajectory beginning 
with an initial offset smaller than the critical value (6j„ < h^ would 
result in a minimum separation smaller than the range of the repulsive 
interaction. Since such separations are unattainable in the presence 
of the interaction, once the particle attains e, the separation remains 
constant and equal to e on the approaching side and, as shown in figure 
|2| the outgoing half of the trajectory collapses onto the trajectory with 
bout = be- On the other hand, for bin > be, the trajectories, unaffected 
by non-hydrodynamic interactions, are symmetric in the absence of in- 
ertia. The single-valued function relating C,min and 6j„ (or bout) is the 
same as that relating e and be- [221 121] Therefore, we can predict the 
dependence of b^ on the pinching gap for a given e by computing ^min 
as a function of 6j„ for different d/D ratios. In conclusion, the critical 
offset is a macroscopic and experimentally measurable variable, which 
serves as a probe/proxy for the range of the repulsive interaction, a 
microscopic effective variable that is experimentally less accessible.^^ 

3.2. Inertia effects on the critical offset. The fore-aft symmetry 
of the trajectories is broken in the presence of inertia (see figure [3]). 
Therefore, in the presence of non-hydrodynamic interactions, we need 
to define two distinct critical offsets: the critical initial-offset b^^m and 
the critical final-offset bc^out, as shown in figure 13} The critical initial- 
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Figure 4. Various trajectories with bin = 2.0, exhibit- 
ing the effect of pinching in the Stokes regime [Re = 
5 X 10~^ and St = 1.11 x 10~^). The axes are nondi- 
mensionahsed by the radius of the particle. The vertical 
dashed lines at x ~ ±7 represent the initial and final 
x-coordinates of the particle, respectively. 

offset bc^in is defined in the same manner as in case of negligible inertia. 
The definition of the critical final-offset (bc^out) needs the auxilliary as- 
sumption that inertia is negligible when the particle approaches small 
surface-to-surface separations from the obstacle. Under this assump- 
tion, the final offset would also present a critical behavior such that, 
all trajectories with bin < bc,in would collapse onto a single trajectory 
downstream to the obstacle, which in turn defines the critical final- 
offset bc,out corresponding to bc^in- 

4. Results and implications of the hard-core model 

4.1. Effect of pinching in the Stokes regime. In figure [4], we show 
representative trajectories corresponding to different apertures of the 
pinching gap, for the same initial offset bin = 2a and with negligible 
fluid and particle inertia. As expected, the trajectories remain fore- 
aft symmetric in all cases. When the aperture of the pinching gap 
is reduced, we observe that the particle follows trajectories that get 
significantly closer to the obstacle. 

In light of the hard-core model discussed before, let us now con- 
sider the effect of the non-hydrodynamic interactions on the trajecto- 
ries through a pinching gap in the Stokes regime. Consider trajectories 
with the initial offset bin = 2a (see figure 111), and non-hydrodynamic 
interactions with a dimensionless range e = eo. Let the critical offset 
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Figure 5. Various trajectories with 6j„ = 2a, exhibiting 
the effect of inertia in the absence of pinching. The axes 
are nondimensionalised by the radius of the particle. The 
vertical dashed lines at x ~ ±7a represent the initial and 
final ^-coordinates of the particle, respectively. 



corresponding to eo in the absence of pinching be be = 6co < 2a. Using 
the arguments in the previous section, in the absence of pinching the 
particle trajectories under consideration are fore-aft symmetric, attain 
a minimum separation C,min > ^o, and the non-hydrodynamic interac- 
tions have no effect on them (e.g., the top-most trajectory in figure [2]). 
As the aperture of the pinching gap is reduced, while keeping the initial 
offset constant at 2a, figure |4] shows that the particle comes closer to 
the obstacle. Eventually, at a particular aperture, we get ^rnin = ^o, 
which implies that the critical offset has increased from bco to 2a. In 
other words, the critical offset that leads to a given minimum sepa- 
ration increases as the aperture of the pinching gap decreases. This 
result is in qualitative agreement with our experiments. [15] In the sup- 
plementary information, [23] we arrive at the same inference directly by 
considering the minimum separation ^min as a function of the initial 
offset bin (see figures SI and S2, and the corresponding discussion in 
the supplementary information |23 j ) . 

4.2. Effect of inertia in the absence of pinching. Figure [5] shows 
representative trajectories (with bm = 2a) in the presence of inertia 
effects, but in the absence of a pinching wall. We observe that the 
effect of inertia is similar to that of pinching discussed earlier, such 
that the particle gets closer to the obstacle as the magnitude of (fluid 
and/or particle) inertia increases. Although the effects of pinching and 
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Figure 6. The dimensionless surface-to-surface separa- 
tion between the particle and the obstacle, as a function 
of the dimensionless x-coordinate of the particle. The 
plot shows the combined effect of particle and fluid in- 
ertia on the dimensionless separation along particle tra- 
jectories for bin = '^CL, in the absence of pinching. The 
vertical dashed line depicts the symmetry plane of the 
system at x = 0. 

inertia on the minimum separation are similar, in the presence of inertia 
the minimum separation is no longer attained at the symmetry plane 
{x = 0) but upstream to the obstacle-center (at negative x values, see 
figure [6]) . We also observe that the x-coordinate at which the minimum 
separation is attained decreases with increasing inertia (similar plots 
showing the individual effects of fluid and particle inertia can be found 
in the supplementary information |23j). 

4.2.1. Effect of pinching and inertia. In figure [7| we consolidate our 
results and present the minimum separation ^^m for all the cases of 
pinching and inertia considered in this work, at bin = 2a. As discussed 
so far, we observe that reducing the aperture of the pinching gap causes 
the particle to attain a smaller minimum surface-to-surface separation 
in all the cases shown, with or without particle and/or fluid inertia. 
Further, the figure also shows that for any fixed aperture of the pinch- 
ing gap the minimum separation attained by the particle decreases 
with increasing inertia. Our simulations also indicate that the effect 
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of particle inertia on ^^m is attenuated in the presence of the pinching 
wall. The attenuation could be due to the hydrodynamic resistance of- 
fered by the pinching wall, which reduces the particle velocity, thereby 
diminishing the effective magnitude of particle inertia (note that for 
calculating the Stokes numbers reported in this work, we have used 
the bulk settling velocity of spheres). The detailed plots of ^min as a 
function of bin for different Stokes and Reynolds numbers are presented 
in figures S2 and S3 in the supplementary information. [23] The effect of 
pinching on these plots is shown in figures S5-S7 in the supplementary 
information. [23] 

We are also interested in the effect of inertia in combination with 
pinching on the final offset b^ut, since changes in the final offset as 
a function of inertia could lead to separation based, for example, on 
density differences. In figure |8} we present the effects of particle and 
fluid inertia on the final offset, for bin = '^cl. In terms of particle 
inertia, we observe that the final offset decreases with increasing Stokes 
number. The final offset decreases less when the trajectories undergo 
pinching, than in the absence of pinching. As mentioned earlier, we 
suspect that this behavior is to the diminished particle inertia due to 
the hydrodynamic resistance offered by the wall. In the case of fluid 
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Figure 8. The final offset bout (corresponding to b^ = 
2a as a function of Stokes (top) and Reynolds (bottom) 
numbers, for different apertures of the pinching gap. 



inertia, we observe that for Re < 0.1 the final offset is unaffected, 
irrespective of the presence or absence of the pinching wall. For higher 
magnitudes of fluid inertia {Re ~ 0(1)), we observe that the final 
offset decreases with the aperture of the pinching gap. Such decrease is 
comparable to the effect of lift forces in the proximity of the pinching 
wall, which render the final offset ill-defined. In the supplementary 
information [23] (figures S5, S6 and the accompanying discussion) we 
present the minimum separation C,min as a function of the initial offset 
bin, which is consistent with the attenuation of particle inertia and the 
presence of lift forces. 

Figure [7] shows that the minimum separation decreases for smaller 
apertures of the pinching gap for a constant magnitude of inertia. 



Therefore, similar to the discussion in ^4.1, we can infer that the crit- 
ical initial-offset bc^in increases with pinching. Further, the minimum 
separation also decreases with an increasing magnitude of inertia, for 
a fixed aperture of the pinching gap. Thus, the critical initial-offset 
increases with increasing particle inertia (see figures S2 and S5-S7 in 
the supplementary information [23]). Since bc,out is a monotonically in- 
creasing function of bc^in, one might expect that bc^out also increases 
with pinching and inertia. However, figure [8] shows that for a constant 
bin, the final offset decreases due to the effect of inertia. Due to such 
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Figure 9. The critical final-offset bc^out corresponding 
to e = 0.3, as a function of Reynolds and Stokes numbers, 
in the presence and the absence of a pinching wall. 



competing effects, commenting on the behavior of the final critical off- 
set as a function of the aperture of the pinching gap and inertia is not 
trivial. Our simulations indicate, as shown in figure [9} that the critical 
final-offset bc,out indeed increases with inertia, but the change becomes 
substantially smaller in the presence of pinching, provided that the ef- 
fect of lift forces generated due to the pinching wall are minimal (also 
see ^min as a function of bout in figures Sll and S12 in the supplementary 
information [23j ) . 

Thus, in conclusion, both critical offsets increase with inertia and/or 
pinching, for a constant range of the non-hydrodynamic interactions 
^min = c- Thus, in a mixture of particles of the same size, a particle 
with sufficiently higher density would exhibit larger critical final-offset 
bc,out than a particle with lower density, thereby rendering its sepa- 
ration possible (assuming that both particles exhibit similar range of 
non-hydrodynamic interactions). In the presence of a pinching wall, 
a similar inference can be drawn. Although, the change in the criti- 
cal offsets due to inertia can be expected to be small because of the 
attenuation of inertia effects by the wall. 
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5. Summary 

We have shown that in the presence of inertia effects and/or a pinch- 
ing wall, suspended particles driven around an obstacle by a constant 
force, get closer to the obstacle due to pinching as well as inertia. The 
minimum surface-to-surface separation C,min as a function of the rel- 
ative pinching d/D quantifies this phenomenon. Augmented with a 
simple hard-core repulsion model for non-hydrodynamic interactions 
(such as solid-solid contact due to surface roughness) our observations 
lead to the definitions of critical initial- and final-offsets. Both offsets 
are equal when inertia is negligible. The offsets increase upon decreas- 
ing the aperture of the pinching gap, and/or upon increasing inertia. 
When the inertia and pinching effects are combined together, the final 
offset varies only weakly with inertia. In this case, we surmise that 
the presence of a pinching wall attenuates the effect of particle iner- 
tia. In conclusion, particles with different densities could be separated 
in systems with significant inertia, more effectively in the absence of 
pinching than in its presence. 
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Susp3d - available to us. This work is partially supported by the Na- 
tional Science Foundation Grant Nos. CBET-0731032, CMMI-0748094 
and CBET-0954840. This work used the resources of the National En- 
ergy Research Scientific Computing Center, which is supported by the 
Office of Science of the U.S. Department of Energy under Contract No. 
DE-AC02-05CH11231. 
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Appendix A. Supplementary material 

A.l. Pinching gap, inertia and minimum separation. Here, we 
present the results of our simulations in a different manner than the 
main text by considering the minimum separation C,min as a function 
of the initial or flnal offsets. We present the results with hierarchically 
increasing the number of variables in the system. First, we show the 
effect of decreasing the aperture of the pinching gap for negligible fluid 



and particle inertia (the Stokes regime, flgure SI). Then we show 



the effect of inertia, in the absence of a pinching wall (flgures S2 and 



S3). Finally we present the cases when fluid and/or particle inertia is 



m 



signiflcant, in the presence of a pinching wall (flgures S5, S6 and S7) 
As discussed in the main text, the relationship between C,Tnin and b. 
shown here, is the same as that between the range of the repulsive non- 
hydrodynamic interactions and the critical offset. [22i i24j Therefore, 
representing the data in terms of ^min^bin plots provides additional 
clarity to the discussion concerning increasing critical offset in the main 
text. 



In flgure S 1 , we present the effect of decreasing the aperture of the 
pinching gap on the minimum separation C,min, in the Stokes regime. 
As noted in the main text, for a given initial offset we observe that ^min 
decreases in the presence of a pinching wall. This is seen by comparing 
the results shown in solid circles with those in open symbols, which 
correspond to the absence and presence of a pinching wall, respec- 
tively. Moreover, comparing the different curves with open symbols, 
which correspond to different apertures of the pinching gap, it is clear 



A SPHERE THROUGH A CONSTRICTION 



19 



0.8 




'1 


• 


- 








• 






- 


....■• 




- 






....• 




- 


0.4 


•■•■■" 




■ O 

. - ■ ■ ■ . -o 


- 


'■ 




..••..□■■ ,.■■■ 


..o ...••'■ 

..A-'' 


- 


0.2 




.■■■■ ..«■■■ 
jo' .■■■ 

.■•■>"' .■■■■■' .•^■' 






e 


o-' ..■■ 


.«■■ ..■•■' 


•■ ■• No wall 
o- o d/D = 0.333 
□■ ■□ d/D = 0.400 
«■ o d/D = 0.500 
A- A d/D = 0.588 




9' 

1 


1 of repulsive forces 


0.1 


1 .- 




, 


- 




1 
1 


1 increases 



1.9 



2.0 

b. 



2.1 



2.2 



Figure SI. The minimum separation C,min as a function 
of the initial offset bin, -Re and St are neghgible. 



0.4 



E 

0.2 



■■°St= 1. 

■ ElSt = 0. 

■■»St = 0. 
■■ASt = 0. 
■■vSt = 0. 
■■•St= 1. 



11x11 

111 

333 

555 

777 

11 



~f Characteristic range 
of repuslive forces 



Initial offset 
increases 



Re = 5x10 



. 




A 


0.8 


V 




.* 






0.4 








0.2 




0.1 







°--ORe = 5xlO 
°--ORe = 5xlO" 
«-»Re = 5xlO 
^■■^Re = 5xl0 



-2 



1.8 1.9 
b. 



St= 1.11x10 



1.6 1.7 1.8 



1.9 
b. 



2.0 2.1 2.2 



Figure S2. The minimum separation ^mj„ as a function 
of the initial offset bin, in the absence of pinching. (Left) 
effect of particle inertia when fluid inertia is negligible, 
(right) effect of fluid inertia when particle inertia is neg- 
ligible. 



20 A SPHERE THROUGH A CONSTRICTION 



0.8 



0.4 



0.2^ 



0.1 



1 


1 


1 ' 1 ' 1 ' 1 

...o 






...•■■■ .■ ...♦ 






..0-- ...•■■ ...■■■ 
..■■■■'■ ..■-»■■■■ 


" 


.0-" 


•■•■....■■•■• ....♦■■■•■ ....■■•* 


- 


..■■■■' 


..♦•■•■' .••*■■ 


o ■ 






■-'' 


..♦■■ 


.>■■■■ 


♦■■■ 








A'' 


o-oRe = 5xlO^St= 1.1 1x10"*, p /P|,= 1 
■■ ■■ Re = 5xl0"^ St = 1. 1 1x10"", p /p^, = 10 
♦■ ■♦ Re = 1.5x10"', St = 1x10"', p /p^. = 30 






A-ARe = 5xl0"',St= l.ll,p /p = 100 


A' 






1 


1 


1,1,1,1 



1.7 1.8 1.9 2 2.1 2.2 

b. 
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that decreasing the aperture leads to smaller minimum separations, 



as the curves shift downwards. Next, in figure [S2] we present the ef- 
fect of inertia in the absence of pinching. The two plots in figure S2 
show the effects of fluid and particle inertia on the minimum separation 
^min, independently. It is clear that ^min decreases when the magni- 
tude of the particle or fluid inertia, increases. The combined effects of 



particle and fluid inertia are presented in figure p3\ where we report 
results corresponding to particle to fluid density ratios 1, 10, 30, and 
100. As explained in the text, inertia causes the particle trajectories 
to be asymmetric. This leads to the minimum separation taking place 



before x = 0. This is corroborated in figure S4 for fluid and particle 
inertia acting independently. It also shows that the x-coordinate at 
which the minimum separation is attained, decreases with increasing 
inertia. 



In figures S5 and S6 , we present the combined effect of increasing Re 
or St for apertures of the pinching gap given hj D = 3d {d/D = 0.333) 
and D = 1.7 d [djD = 0.588), respectively. The plot on the left in 
each figure shows the effect of increasing St for a fixed aperture of the 
pinching gap and Re. The plot on the right in each figure shows the 
effect of increasing Re for the same aperture of the pinching gap and 
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Figure S6. The effect of inertia on the imin-bin rela- 
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constant St. For bm = 2a, the corresponding data is presented in fig- 
ure 7 of the main text. We observe in both cases, that the minimum 
separation decreases with increasing inertia. In the case of a constant 
Re (i.e., increasing St), we see that decreasing aperture of the pinching 
gap (comparing figures S5 and S6) attenuates the effect of inertia on 
the minimum separation for a fixed initial offset. The decrease in mini- 
mum separation (between two magnitudes of inertia) for a given initial 
offset is smaller when the aperture of the pinching gap is decreased. In 
contrast, the effect of increasing Re seems to increase as the aperture 
of the pinching gap becomes smaller, specifically for Re ~ 0(1). In 
this case, however, we observe significant lift due to particle-wall in- 
teraction, which makes the discussion in terms of offsets less relevant. 
Note that for Re ^ there are no asymptotic offsets, and we observe 
that for Re ~ 0(1) the effect of lift becomes significant over the length 
of the system. 

Finally, representative examples of the combined effect of varying all 
Re, St and d/D are shown in figure S7 The plot on the left shows 



the results corresponding to four pairs of St and Re for d/D = 0.333. 
These pairs correspond to the density ratios mentioned earlier, Pp/pf = 
1, 10, 30, 100, for a constant force acting on the particle. The plot on 
the right presents the minimum separation for the same pairs of Re 
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Figure S7. Combined effect of tlie pinching wall, par- 
ticle and ffuid inertia on tlie ^min^bin relationsliip. (Left) 
tlie aperture of tlie pincliing gap is 2)d (i.e., d/D = 0.333), 
(right) the aperture is 1.7d (i.e., d/D = 0.588). 

and St but for, d/D = 0.588. The results follo"w that, increasing iner- 
tia effects (both particle and fluid inertia) leads to smaller minimum 
separations for a flxed initial offset. Also analogous to the results dis- 
cussed earlier, increasing the extent of pinching decreases the minimum 
separation corresponding to a given initial offset, but attenuates the 
relative reduction in minimum separation due to inertia effects. 



A. 2. Minimum separation, critical offsets and non-hydrodynamic 
interactions. Similar to the discussion in the main text, here we 
present the effect that inertia and a pinching wall have on the criti- 
cal flnal-offset of the particle trajectories. It is clear from the results 
discussed both in the main text as well as in the previous section that 
the minimum separation ^^m decreases when the aperture of the pinch- 
ing gap is reduced, or the magnitude of inertia is increased, or both. 
Using the ^min^bin plots introduced in the previous section, we elabo- 
rate more on the idea of the critical initial-offset bc,in, and the effect of 
pinching and/or inertia on it. Further, extending the method, we plot 
^min as a function of the flnal offset bout for various cases of pinching 
and inertia, and comment on the behavior of the critical flnal-offset 

"cMUt- 



Consider flgure SI in light of the hard-core model of the non-hydrodynamic 
repulsive interactions explained in the main text. For a flxed range of 
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the interactions e, we have explained in the main text that the criti- 
cal final-offset completely characterizes the lateral displacement of the 
particle downstream of the obstacle. Given that the case presented in 
figure SI corresponds to negligible inertia (Stokes regime), both critical 



initial- and final-offsets are equal. Thus, figure SI can also be inter- 
preted as a plot of e versus be, and clearly shows that the critical offset 
be increases with decreasing aperture of the pinching gap. Similar con- 
clusions can be drawn concerning the critical initial-offset in absence of 



a pinching wall using figures [S2] and [S3l that is, the critical initial-offset 
increases for a fixed range of the repulsive interactions, with increasing 
magnitude of inertia (fluid and/or particle inertia). 

The asymmetry imparted to the particle trajectories by particle iner- 



tia, is readily seen in figure S8 This figure supplements the information 
provided in figures 3 and 5 in the main text. We observe that, in the 
absence of a pinching wall, the final offset decreases with increasing 
particle inertia, while fluid inertia has a weak influence on it. As a 
result, the effect of inertia competes with the effect of pinching and in- 
ertia (decreasing the final offset for a given initial offset), as discussed 
in the main text (in the context of the dependence of bc^out on pinching 
and inertia). Figures S9 and SIO present the effect of inertia on the 
minimum separation, as a function of the final offset, and complement 
figures S2 and S3 respectively. These figures show that the critical 
final-offset increases with the magnitude of inertia for a constant range 
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Figure S9. The ininimum separation as a function of 
tlie final offset, (left) effect of particle inertia, (right) ef- 
fect of fluid inertia. 
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Figure SIO. Effect of both significant fluid as well as 
particle inertia on the ^min^bout relationship, in the ab- 
sence of pinching. 



of the repulsive interactions, e. The presence of both pinching as well 
as inertia poses the most complex case for analysis, especially, when 



non-hydrodynamic repulsive forces also play a role. Figures SIX and 
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Figure S 1 1 . The individual effect of particle and fluid 
inertia in the presence of pinching, on the ^min-bout re- 
lationship. The aperture of the pinching gap is 3.0(i 
(i.e., d/ D = 0.333). (Left) The effect of particle iner- 
tia, (right) the effect of fluid inertia. 
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Figure S12. The individual effect of particle and fluid 
inertia in the presence of pinching, on the imin^bout re- 
lationship. The aperture of the pinching gap is 1.7d 
(i.e., d/D = 0.588). (Left) The effect of particle iner- 
tia, (right) the effect of fluid inertia. 
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Figure S13. The effect of inertia on the C,min-bout re- 
lationship, in the presence of pinching for two different 
apertures of the pinching gap. (Left) The aperture of 
the pinching gap is S.Od (i.e., d/D = 0.333), (right) the 
aperture is 1.7d (i.e., d/D = 0.588). 



S12 are analogous to figures S5 and S6 Again, we see that the effect 



of pinching is to reduce the minimum separation and that the effect of 
inertia is small, except for Re ~ 0(1), where we observe a reversal of 
the usual trend, with high inertia leading to smaller values of the final 
offset. As discussed in the main text, this is caused by lift forces, that 
reduce the final offset, and is therefore a case that should be considered 
separately. 



